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Abstract 

We study ghosts in multimetric gravity by combining the mini-superspace and the Hamiltonian 
constraint analysis. We first revisit bimetric gravity and explain why it is ghost-free. Then, we 
apply our method to trimetric gravity and clarify when the model contains a ghost. More precisely, 
we prove trimetric gravity generically contains a ghost. However, if we cut the interaction of a pair 
of metrics, trimetric gravity becomes ghost-free. We further extend the Hamiltonian analysis to 
general multimetric gravity and calculate the number of ghosts in various models. Thus, we find 
multimetric gravity with loop type interactions never becomes ghost-free. 
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I. INTRODUCTION 



It is interesting to explore the possibility that a graviton is massive both from theoretical 
and phenomenological point of view. Theoretically, it is challenging because of various 
difficulty in constructing a consistent model for a massive graviton. At the linear level, 
Fiertz and Pauli succeeded in constructing a ghost-free model for a massive graviton 
However, it is soon recognized that there is a tension between the theory and experiments, 
the so-called van Dam-Veltman-Zaharov (vDVZ) discontinuity 

a a. 

It is suggested that 

the non-linearity resolves the vDVZ discontinuity [4!]. Unfortunately, it turned out that the 
non-linearity gives rise to a ghost, the so called Boulware-Deser (BD) ghost sl. Recently, 
de Rham, Gabadadze and Tolly have succeeded in constructing ghost-free non-linear massiv 
gravity theory {g, 7| (See a review [8] and references therein.). Still, there remains various 
theoretically intriguing issues to be explored. Phenomenologically, there is a chance to 
explain the current accelerating universe based on massive gravity. In fact, there appears 
an effective cosmological constant proportional to the square of graviton mass UH]. It IS 
worth studying this possibility in detail. 

One peculiar feature of massive theory of graviton is the necessity of a reference metric 
which breaks the diffeomorphism invariance. It is natural to promote this reference metric 
to a dynamical variable, which is nothing but bimetric gravity. Bimetric gravity contains 



two metrics g and / interacting each other. A history of bimetric gravity is long |15Nl7|. 
Curiously, bimetric gravity also suffers from the ghost problem. Thanks to the recent devel- 
opment in massive gravity, however, Hassan and Rosen have proposed ghost-free bimetric 



gravity 18l-l20|. Then, a natural question arises whether or not we can construct ghost-free 
multimetric grav ity. Actually, a naive extension of bimetric gravity to trimetric case was 
proposed in 21']. There, three metrics f and h have a pair interaction between {g,h) 
,{h, f) and {f,g), which forms a loop structure. In contrast to bimetric gravity, however, 
the presence or the absence of BD-ghost remains unknown. Recently, Hinterbichler and 
Rosen showed that a large class of multivielbein gravity is ghost free 22|. The relation to 
metric theory is also discussed |23|. However, the relation to the models presented in 21 1 
is not clear. The difference between 21 1 and {22! comes from the loop type interaction. In 
fact, the proof by vielbein method is not applicable to the loop type interaction. Hence, we 
need to study multimetric gravity with a different approach. 



In this paper, we propose a simple method to study the ghost problem and clarify when 
multimetric gravity is ghost-free. A method often used for the ghost analysis is to examine 
models in the decoupling limit. However, a more honest way for probing ghosts is to use the 



Hamiltonian constraint analysis using the ADM formalism [2^. The difficulty in studying 
multimetric gravity with the constraint analysis comes from the existence of a shift vector. 
To avoid the difficulty, we employ the mini-superspace approximation. The mini-superspace 
reduction of phase space makes the analysis so simple. Nevertheless, it is sufficient to identify 
ghosts because this reduction process does not fail to capture ghosts. 

The organazation of the paper is as follows. In section [Tll we revisit bimetric gravity and 
explain our strategy for the ghost analysis. In section IIIIl we investigate trimetric gravity 
using our method and found that the loop type interaction allows ghost. In section [TV] we 
further extend the analysis to general A/'-metric gravity. We clarify when ghost appears in 
the spectrum. The final section |V] is devoted to the conclusion. 



II. BIMETRIC GRAVITY REVISITED 

In this section, we revisit bigravity and explain our method to probe a ghost. It is already 
known that bimetric gravity is ghost-free [l9|. Here, we show the same conclusion can be 
obtained using a simple mini-superspace approximation. In the context of massive gravity, 
the decoupling limit analysis turns out to be a useful way for the ghost analysis. However, 
the most complete one is to use Hamiltonian constraint analysis and count physical degrees 
of freedom. Our strategy is to use Hamiltonian constraint analysis in the mini-superspace. 

The action of ghost-free bimetric gravity {20 1 is given by 

= Ml j d^x^-deigR[g\ + M] j d^x^-dei fR[f] 

r ^ 

+2mX/ / d^x^-deigY,(^nen{^/9^f), (1) 

n=0 

where the first and the second terms are Einstein-Hilbert action for each metric g and / 
from which we can calculate the scalar curvatures R[g]i R[f]- Here, we have two Planck 
masses Mg and Mj. The last term describes the interaction between two metrics and /3„ are 
dimensionless coupling constants. The other constants m and Mgj are introduced to adjust 
the mass dimension. The square root of the matrix is defined such that \/ g^^j \/ 9^^J = 



(D KD 



FIG. 1. The interaction \/g can be represented by a simple diagram. Each blob describes a 
spacetime with a given metric. The arrow indicates the order of product. 

g^^fxv The interaction terms are constructed by e„(X) which we define, for matrix X, 



eo(X 
ei(X 
62 (X 
e3(X 

64 (X 



trX 

^ (tr^X - trX^) (2) 
i (tr^X - StrXtrX^ + 2trX3) 

— (tr^X - etr^XtrX^ + Str^X^ + StrXtrX^ - 6trX^) = detX , 
24^ ' ' 



where we used the notation tr"X = (trX)" and trX" = tr(X"). It is useful to represent the 
interaction by a diagram in Fig. [TJ Note that there is the order between g and / which 
is denoted by the arrow. It is known that the interaction produces massless and massive 
gravitons and the spectrum is free of Boulware-Deser ghost. This feature comes from a 
specific interaction form found in massive gravity theory. Remarkably, there exists the 
diagonal diffeomorphism invariance which indicates the presence of the massless graviton. 

Now, let us perform Hamiltonian constraint analysis based on the ADM formalism. In 
particular, to make the analysis tractable, we employ the mini-superspace approach. Namely, 
we assume spatial homogeneity and express metrics in terms of ADM variables as 

g^udx^'dx'' = -Nitfdt^ + ^ij{t)dx'dx^ , (3) 

where is a lapse function and is a spatial metric. Similarly, we can take the following 

ansatz 

f^^dx^'dx" = -L{tfdt^ + Uij{t)dx'dx^ , (4) 
where L is a lapse function and Uij is a spatial metric. It is convenient to write them in a 
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matrix form, 





u^=\ : " I , r=\ ^'^ I , (5) 



where 7'-' and w'-^ are inverse matrices of spatial metrics and . Then, a basic part of 
interaction terms can be calculated to be 



i'ui, ] V 



7 ^ijj 



When we count physical degrees of freedom, the following must be taken into account. 
In the vacuum cases, we can diagonalize one of two spatial metrics using diagonal spatial 
coordinate transformations. Performing a spatial coordinate transformation — > A(to)*j 
we can set one spatial metric at the time t = to, 7y(to), a unit matrix 6ij. Moreover, since 
the orthogonal transformation dose not change 7ij(to) = ^ij, we can diagonalize 7jj(to) 
simultaneously by using this freedom. At this stage, homogeneous spatial coordinates are 
completely fixed. Now, jij and 7^ is diagonal at the time t = to as an initial condition. 
Then we assume diagonal form of 7ij(t) at all time, and insert it into equations obtained 
from variations of action. Any contradiction never occurs in vacuum. Thus, we conclude 
that one spacial metric 7ij(t) can be diagonalized because of the uniqueness of the solution. 
Hence, the number of component of one of two metrics reduces from 6 to 3. This fact will 
be used later. 

In this paper, for simplicity, we assume that interactions are minimal 
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25| . namely 



/3o = 3, A = -1, /32 = 0, /33 = 0, /34 = 1 . (7) 

Clearly, this simplification does not lose any generality concerning with the ghost analysis. 
Then, the Lagrangian reads 

£ = My^'y,^ + M}p'^uj,, - NCn - LCl , (8) 
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where tt*-' , p''^ are canonical conjugate momentum of jij, ouij . Here, we have defined 

M2 



Um — 



■■ {^''^iJ - - K ^'^^[7] + ai Vd^(tr V^P^ - 3) (9) 



-\/det7 

Cl = - ]:P\p'^ - M?^d^ (3)i?[cc;] + ai (^d^ - Vd^) , (10) 

Vdetu;^ 2 V 

where the first two terms of constraints and Cm come from Einstein-Hilbert term in 

the action, so *^^)i?[7] and ^^^i?[a;] are spatial scalar curvatures computed from 7 and a;, 

respectively. The last term of each constraint comes from the interaction (see the derivation 
in Appendix A), and we use oi = 2w?Mg^. Now, the Hamiltonian is given by 

H = NCn + LCl . (11) 

Since there are two Lagrange multipliers, there are two primary constraints 

Cat = 0, = . (12) 

Moreover, we need to impose consistency conditions for them 

Cn — {Cn , i?} — L{Cn , Cl} = LCnl ~ , 

Cl = {Cl,H} = N{CL,Cr,} = NCln ~ , (13) 

where the Poisson bracket {F,G} is defined by 

{F C] -( \ ( \ (U) 



Here, "~ 0" means "= 0" on the constraint surface. Notice that {F,F} = because of 
spatial homogeneity. 

To check if secondary constraint arises or not, we have to calculate Poisson bracket C^l- 
Prom the calculation presented in Appendix B, we obtain 



Cnl — {Cn , Cl} — 



ai 



2 My, - ^ ^ ( tr - tr ( p a;) 



• (15) 



This leads to one secondary constraint Catl ~ 0. The consistency condition for the secondary 
constraint reads 

Cnl = N{Cnl , Cn} + L{Cnl ,Cj}^0. (16) 



This condition determines one of two Lagrange multipliers N and L. The remaining multi- 
plier describes the diagonal time reparametrization invariance in bimetric gravity. 

The number of components of two metrics and their canonical conjugates is 24. Since 
we can diagonalize one of the two metrics, we should subtract 6 from this number. Recall 
that there are two primary constraints and one secondary constraint. Furthermore, as we 
have one first class constraint, we have to put one gauge condition. Thus, the total number 
of degrees of freedom should be (24 — 6 — 2 — 1 — l)/2 = 7in configuration space, which 
matches degrees of freedom of one massless graviton and one massive graviton. This proves 
that BD ghost is absent in bimetric gravity described by the action (|T]). 



III. TRIMETRIC GRAVITY 



Now, we apply the method explained in the previous section to trimetric gravity. In 
contrast to the bimetric gravity, there are two kind of interactions, namely, the tree type 
and the loop type interactions. We discuss both cases, separately. 

The action for trimetric gravity [21] can be written as 

Stri = Ml j d^x^- det gR[g\ + M] J d^x^ - det fR\f\ + J d'^xV- det hR[h] 

r ^ 

+2m?M|^ / d^x ^/-detg ^ /3„e„ ( a/^) 

n.=0 
n=0 

+2mlMl / d^xV^d^J2^'M\^)^ (17) 

n=0 

where (3'^ and are free parameters and R[g], R[f] and R[h] are scalar curvatures 
constructed from metrics g, f and h, respectively. We also introduced new mass parameters 
7711,1712,171^, Mfh, Mhg and a Planck mass Mh- It should be noted that there exists the 
diagonal diffeomorphism invariance in this trimetric theory which makes one of gravitons 
massless. As discussed in 



211 . |22|, if this trimetric gravity contains no extra degrees of 
freedom, the total number of degrees of freedom should be 2 + 5 + 5 = 12, which comes from 
one massless graviton and two massive gravitons. From now on, we use 

a, = 2mlMlf, a2 = 2mlM%, = 2mlMlg (18) 
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for notational simplicity. If we have oi 7^ 0, 02 7^ and 03 7^ 0, all pairs {g, /), (/, h) and 
{h,g) interact and we call it the loop type interaction. When one of {i = 1,2,3) is set 
to zero, two of three pairs of interactions remain, which we call the tree type interaction. 
The case where one interaction is cut is already proved to be ghost-free using vielbein 
formalism 22|, however, for the loop type interaction no one shows the presence or the 
absence of ghost. In this paper, we settle this issue. 

Apparently, the full Hamiltonian constraint analysis is difficult. To circumvent this dif- 
ficulty, we take the method used in the previous section. Namely, we assume spatial homo- 
geneity and express metrics in terms of ADM variables as 

gf.^dx^dx" = -Nitfdt^ + -fij{t)dx'dx^ , (19) 

where iV is a lapse function and 7ij is a spatial metric. Similarly, we can take the following 
ansatz 

ff.^dx'^dx" = -L{tfdt^ + Uij{t)dx'dx^ , (20) 

and 

h^^dx^'dx" = -Q{tfdt^ + pij{t)dx'dx^ , (21) 

where L and Q are lapse functions and Uij and pij are spatial metrics. To perform Hamil- 
tonian constraint analysis, we need the Lagrangian in the ADM variables 

C = My^% + M}p'^uj,, + Ml<P'^p,^ - NCm - LCl - QCq , (22) 

where vr'-', p^^ and 0*-^ are canonical conjugate momentum of 'jij, Uij and pij. Here, three 
Hamiltonian constraints 

+aiA/det7(trA/7~iw - 3) + a3{^/det p - ^/det'y) , (23) 



V det u ^ ^ ^ 



0J\ 



and 



+a2V^deta;(trA/a;^^p — 3) + ai(A/det7 — Vdetw) , (24) 



-t-as^/det p{i'c\J p^^^ — 3) + 02 [V det u — \/ det p) (25) 
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emerge. The first line of each Hamiltonian constraint comes from the Einstein-Hilbert term 
in the action, so ^^^i?[7], and are spatial scalar curvatures calculated from 7, 

uj and p, respectively. The other terms can be derived as explained in Appendix A. Then, 
the Hamiltonian can be read off as 

H = NCn + LCl + QCq . (26) 

Since there are three Lagrange multipliers, there arise three primary constraints 

Civ = 0, Cl = 0, Cq = Q. (27) 

Moreover, we need consistency conditions for them 

Cn = {Cat , H} = L{Cn , Cl] + Q{Cn , Cq} f« , 

Cl = {Cl , H] = N{Cl , Cm] + Q{Cl , Cq} ^ , (28) 

Cq = {Cq , H] = N{Cq , Cm] + L{Cq .Cl]^Q. 

To check if secondary constraints arise, we must calculate Poisson brackets. From the 
calculation in Appendix B, we obtain 



Cnl = {Cn , Cl] — ai 



'■'^ V deta; V2^ 

By performing permutations among g = {N, 7) , / = {L, u) and h — {Q, p), we also get 

^1 



• (29) 



ClQ — {Cl , Cq] — 02 



' det UJ /I 



det p 



trVc^-V-tr(Vc^-V0p)) 



(30) 



and 



Cqn = {Cq , Cn] — (I3 



.(31) 



In general, quantities inside the bracket does not vanish. Hence, the coefficients 01,02 and 
03 determine the consistency conditions. 

A. Tree Type Interaction 

In this subsection, we consider the tree type interaction 



oi 7^ 0, 02 = 0, 03 7^ , 



(32) 



KD 



FIG. 2. The diagram represents the tree type interaction. 

which cut interaction between / and h as in Figj2l In any case, there are primary constraints 
fl27|) . Since Cql = Clq = trivially holds, consistency conditions f l28|) lead to equations 

LCnl + QCnq^O, NCln-0, NCqn^O. (33) 

Hence, we have two secondary constraints 

Cnl ~ 0, Cnq^O. (34) 

Moreover, we must impose consistency conditions 

Cnl = {Cnl , = N{Cnl , Cjy] + L[CiyL , Cl] + Q{Cnl , Cq] ^ 

Cql = {Cql,H}= N{Cql , Cn] + L{Cql , Cl} + Q{Cql ,Cq}^0, (35) 

which determine two of three Lagrange multipliers A^, L and Q. The remaining multiplier 
is related to the gauge transformation. 

Eventually, we have five constraints and one gauge freedom. In trimetric gravity, prop- 
agating modes are spatial metrics. Each of them has six components, but as is already 
explained we can diagonalize one of them. Hence, trimetric gravity has 3 + 6 + 6 = 15 
degrees of freedom in configuration space and 15 x 2 = 30 in phase apace. Thus, the total 
number of degrees of freedom is (30 — 5 — 1)/2 = 12 which matches the physical degrees 
of one massless and two massive gravitons. Therefore, no BD ghost exists in the spectrum. 
This conclusion is consistent with the one obtained by the vielbein method [2^. 
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FIG. 3. The diagram represents the loop type interaction. 



B. Loop Type Interaction 

Now, we consider the more general loop type interaction represented by a diagram in 
Figl 

ai ^ 0, as ^ 0, as ^ . (36) 

It is obvious that 

[Cn.Cl]^^, {Ci,CQ}^0, {Cq,C^}^0 (37) 

even on the constraint surface. Hence, consistency conditions ( l28l) do not generate any 
secondary constraint. Instead, it determines Lagrange multipliers A^, L and Q. However, 
due to the antisymmetric property of Poisson brackets 

Cnl = {Cn , Cl} = —{Cl , Cn} = Cln (38) 
Clq = {Cl , Cq} = -Cql (39) 
Cqn = {Cq , Cn} = —Cnq, (40) 

only two of them are determined. For example, choosing 

L = -^Q, N = -^Q, (41) 

all of consistency conditions ( |28|) are satisfied. 

To conclude, we have three primary constraints and we need one gauge condition to fix 
one undetermined Lagrange multiplier which is associated with the time reparametrization 
invariance. In trimetric gravity, as is already counted, there are 3 + 6 + 6 = 15 degrees 
of freedom in configuration space and 15 x 2 = 30 in phase apace. In phase space, we 
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have three constraints and one gauge condition, so total number of degrees of freedom is 
(30 — 3 — l)/2 = 13. If no BD ghost is present, there must be 2 + 5 + 5 = 12 degrees 
of freedom which comes from one massless graviton and two massive gravitons. Therefore, 
one extra degree of freedom exists and it should be a BD ghost. Thus, we have proved the 
existence of a ghost in generic trimetric gravity. 

IV. GENERAL MULTIMETRIC MODELS 

Now, we are in a position to discuss more general cases. We exphcitly calculate the 
number of ghosts if they exist. 

In this section, we consider M dynamical metrics Qk {k = 1, 2, .., TV) and interaction terms 
such as 

Af 4 

ak V- det Qk ^ h,nen(\J Qk^Qk+i) , (42) 

k=l n=0 

where we define 5'a/-+i = gi and for later purpose we also need go = g^. Let us describe the 
interaction between two metrics g^ and gk+i in terms of ADM form of metrics 

dsl = -Nl{t)dt^ + -fk,ij(t)dx'dx^ . (43) 

Schematically, the interaction can be written as 

4 

V- det gk ^ /3k,ne„ (^y gk^gk+ij = NkFk{jk ■ 7fe+i) + Nk+iGk{jk ■ Ik+i), (44) 

n=0 

where Fk and Gk are some functions determined by parameters f3k,n- Thus, the total inter- 
action terms are given by 

M 4 M 

Qky/- det gk f3k,nen (\j g^^gk+i] = A^fe{afe-^fe(7fe : 7fc+i) +0^-1^^-1(7^-1 ^ 7fe)} ■ 

fe=l n=0 fe=l 

(45) 

The Hamiltonian becomes 

H ^Y NkCk, Ck = C|g(7fe, TTfe) - akFk{-ik ■ 7fe+i) - afe-iGfe-i(7fe-i : Ik) , (46) 

k=l 

where comes from the Einstein Hilbert term for gk, so it contains ^k and its canonical 
CO jugate momentum tt^. 
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Corresponding to J\f Lagrange multipliers, we have M primary constraints 



Cfc = 0, (A; = 1,2,..,AA) . 



Next, we have to examine N consistency conditions 



Ck — Ck,k-iNk-i + Ck,k+iNk+i ~ 0, 



(47) 



(48) 



where Ck,i = {Ck, Q} and Ck,i = if |/c — /| ^2. In this formula, A^o = and 
Nj^j^i = Ni should be understood. Note that the explicit calculation gives rise to an 
important information 



Ck,k+i oc ttk . 



(49) 



The structure of this matrix depends on odd or even number. For example, in the case 
jV = 4, we have 



Ck,i — 



( Ci,2 Cia\ 

— Ci^2 (^2,3 

-C2,3 C3,4 

V-Ci,4 -C3,4 / 



(50) 



While, in the case of = 5, we get 

/ Ci,2 Ci,5\ 

-Cl,2 C2,3 

Ck,l= -C2,3 C3_4 . (51) 

-C3,4 Q,5 

V-Ci,5 -C4,5 I 

In the case of odd number of metrics, we cannot split the equations into two independent 
sets. While, in the case of even number of metrics, we can split a set of equations into 
independent two groups of equations. Hence, we have to discuss two cases, separately. 



A. Tree Type Interaction 

First, we consider the tree type interaction. 
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FIG. 4. The diagram represents the tree type interaction. 

If we cut one of {gk, gk+i)k=i,2,..,M interactions as in FigJU for example setting ai = 0, 
Eq. ( I48l) leads to A/" — 1 secondary constraints 



Ck,k+i^O, (/c = 2,3,..,A^), 



(52) 



and their consistency conditions 



Ck,k+i = J2 i^^-'^+i ' ^'}^' ~ 0, (A; = 2, 3, .., AT) 
1=1 



(53) 



determine A^ — 1 of A^^ {k = 1,2, .., A/"), only one Lagrange multiplier remains undetermined. 
Therefore, the total number of degrees of freedom can be deduced as 



^(2(3 + 6(Ar-l)) -A/--(A/--i)-i) =5(Ar-l) + 2, 



(54) 



which corresponds to A/" — 1 massive and one massless gravitons. Therefore, there exists 
no BD ghost. This conclusion is also consistent with the one obtained by the vielbein 



method 



22|. 



B. Loop Type Interaction 

Now, we come to our main point. 

If all of {gk, gk+i)k=i,2,..,Ar interactions exist as in FigJSl the analysis gets a little compli- 
cated. We have to discuss odd and even numbers, separately. 
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FIG. 5. The diagram represents the loop type interaction. 



1 . Odd Number of Metrics 



First, we consider the case where M = 2m + 1, where m is a natural number. In this 
case, we can classify Eq. fHHj) into the following four parts 



C2k,2k-lN2k-l + C2k,2k+lN2k+l — 
C2fc-l,2fc-2^2fc-2 + C2k-l,2kN2k = 

,2m+1^2m+l + (^1,2^2 — 
C2m+l,2m^2m + C2m+l,l^l = . 



(A: = l,2,3,...,m) (55) 
(A; = 2,3,...,m) (56) 

(57) 
(58) 



Solving Eq. (j55ll . we see all of A^2fc+i {k = 1,2, ..,m) can be expressed by Ni. Similarly, 
Eq.( l56i) can be used to express N2k {k = 2, 3, .., m) in terms of N2. Substituting these results 
into Eq. (|5711 and Eq.( |58l) . we obtain a single equation which determines N2 by A^^i. Thus, 
Eq. (H8!) determines A/" — 1 Lagrange multipliers, and one multiplier is left undetermined, 
which reflects the existence of gauge symmetry. 

In the case of odd number of metrics, there is no secondary constraint. While, we need 
one gauge condition to fix the gauge degree of freedom. In conclusion, the total number of 
degrees of freedom can be calculated as 

^(2(3 + 6(Ar-l)) - AT - 1) =5(A/'-l) + 2 + :^^^^ . (59) 

Here, the first two terms correspond to massive and massless gravitons, respectively. The 
last one should be BD ghosts and the number of ghosts is given by (A/" — l)/2. 
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2. Even Number of Metrics 

Next, we consider the case A/" = 2m + 2, where m is a natural number. In this case, we 
can spht Eq.( H5]) into two independent sets of equations, 

Ck,k-iNk-i + CkMiNk+i = (A; = l,3,5,.,2m + l) (60) 
Ck,k-iNk-i + Ck,k+iNk+i = (A; = 2, 4, 6, .., 2m + 2). (61) 

The first set fl60p contains only A^^ (A; = 2, 4, 6, .., 2m + 2), and the second set fl6ip contains 
{k = l,3,5,..,2m + 1). Here, if the component Ck^k±i is in Eq.f lBU]) . Ck±i,k = —Ck,k±i 
must be in Eq. (16 1 p and vice versa. Therfore, in each set, every component Ck,k±i appears 
only once. Now, we define 

= C2^-l,2j , M, = N2j («,j = l,2,3,..,m+l). (62) 

Note that Dij ^ only for i — j = 0, 1. Then, Eq.( BU]) can be written as 

'^DijMj=0 (i = 1,2,3,. .,m+ 1), (63) 

j 

which we can split into 

Di,iMi + Di,^+iM„+i = (64) 
A,i-iMi_i + A,iMi = (z = 2,3,..,m+l). (65) 

Using Eq. fl65p . we can solve all of Mj {i = 2,3, ..,m + 1) in terms of Mi. However, the 
relation between Mi and M^+i obtained from Eq. (l^ is not the same as Eq. flM|) because 
Eq. (l65!l contains no Di^i and Di^m+i- So, we have to impose a constraint so that we get 
non-trivial Lagrange multipliers. This is a secondary constraint expressed by 

det Dij = . (66) 

Under this condition, m of Mj (j = 1, 2, .., m + 1) are determined, and one is left undeter- 
mined. 

Now, we take latter set (IMj) and define 

Eij = C2i,2j-i , Wj = N2j^i (z,j = 1,2,3,. .,m + l). (67) 
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The same argument applies, so we get a secondary constraint det Ey = 0, and one of 
Wj {j = l,2,..,m + 1) is left undetermined. However, matrix Eij satisfies Eij = —Dji. 
Hence, detEij = is not a new constraint. Therefore, from Eq. fl48l) . we get one secondary 
constraint det Dij = and two undetermined Lagrange multipliers . Then, we must impose 
a consistency condition for the secondary constraint 



which reduces the number of undetermined Lagrange multipliers from two to one. 

To summarize, there are A/" primary constraints and one secondary constraint and we need 
one gauge condition. Thus, we come to the conclusion that the total number of degrees of 
freedom is 



Here, again, the first two terms correspond to massive and massless gravitons, respectively. 
Hence, the number of BD ghosts should be (A^ — 2)/2. 

C. More General Diagrams 

In the previous sections, we have considered tree and loop type interactions. In the 
case of bimetric gravity, the interaction type is unique, namely, there is only the tree tj^e 
interaction. In the case of trimetric gravity, there are two possibilities, the tree and the loop 
type interaction. In the case of tetrametric gravity, there are many loop type interaction 
represented by a diagram (a) in Figj6l If we cut some of the interaction, we can make the 
tree type interaction and the broom type interaction represented by a diagram (b) in Figj6l 
From our analysis, it is apparent that if the interaction contains at least a loop, then there 
are ghosts. For example, the model with a diagram (c) in Figj6] contains a ghost. Therefore, 
in generic cases, there exist ghosts in multimetric gravity. The number of ghosts depends on 
the interaction pattern. To construct a viable model, we have to eliminate all of loop type 
interactions. 




(68) 



k=l 




(69) 
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M (b) tc) 

FIG. 6. The diagram (a) represents the most general type interaction. The diagram (b) is the 
broom type interaction. The diagram (c) includes the loop type interaction, hence there should be 
a BD ghost. 



V. CONCLUSION 



We studied multimetric gravity by combining the mini-superspace and the Hamiltonian 
constraint analysis. We first revisited bimetric gravity and explained why it is ghost-free. 
This proved validity of our method. Then, we applied our method to trimetric gravity and 
clarified when the model contains a ghost. We proved trimetric gravity generically contains 
a ghost. However, if we cut the interaction of a pair of metrics, trimetric gravity turned 
out to be ghost-free. We further extended the Hamiltonian analysis to general multimetric 
gravity and calculate the number of ghosts in various models. Thus, we found multimetric 
gravity with loop type interactions never becomes ghost-free. The number of BD ghost in 
JV metric case turned out to be (A/" — l)/2 or (A/" — 2)/2, depending on whether the number 
of metrics Af is odd or even. Hence, the number of BD ghosts increases by one everytime 
two more metrics are introduced. There are other models which may contain ghosts or may 
not contain any ghost. It depends on the interaction type. The number of ghosts can be 
calculated once the diagram characterizing the interaction pattern is given. 

Admittedly, what we have investigated is BD ghosts. There may be other ghosts de- 
pending on the solutions 0, l^^] . In other words, the absence of BD ghosts is a necessary 
condition as a healthy model. In this paper, we have studied interaction terms consists of 



only pairs of metrics. However, as in 22|, interactions of triplets or quadruplets may be 



allowed. We hope to study this possibility in future. It is also interesting to extend our 



analysis to higher curvature theories 



28|. 



18 



ACKNOWLEDGMENTS 

This work was supported in part by the Japan Society for the Promotion of Science 
(JSPS) grant No. 24 -1693, the Grant-in-Aid for Scientific Research Fund of the Ministry 
of Education, Science and Culture of Japan No. 22540274, the Grant-in-Aid for Scientific 
Research (A) (No.21244033, No.22244030), the Grant-in-Aid for Scientific Research on In- 
novative Area No. 21111006, JSPS under the Japan-Russia Research Cooperative Program, 
the Grant-in-Aid for the Global COE Program "The Next Generation of Physics, Spun from 
Universality and Emergence". 



Appendix A: Interaction Terms 

In this appendix, we calculate the interaction terms separately. We use the following 
representations 

The first one is given by 

^-detgei{y/g-^f) = ^-deigtry/g-^f 

= ydet7 n(^L/N + tr Vt^) 

= Vdet7 (^L + Ntry/-f-^u?j . (A2) 

Due to the combination g^^f, we got the linear terms with respect to the lapse functions. 
The second one becomes 

A/^d^e2(yP7) = v'^d^^(tr2yP7-tr(^-V)) 

= Vdct7 ^ I (^L/N + tr Vt^) ^ - L^N^ - tr (7-^6^) } 

= Vdet7|LtrVo^+ (tr^y7^-tr(7"^cj))| . (A3) 
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Again, we obtained desired linearity for the lapse functions. The third one can be calculated 

as 



V^^d^ i (tr^ - StTy^tT{g-'f) + 2tT{g-'ff' 

~'i[L/N + trv^^^) [l^/N^ + tr(7-ia;)) + 2(^LViV=^ + tr (7-^0;)'/') 
v/det7 I (ti^ v^T^ - tr (7-^0;) ) 

+ ^A^(^tr3v/^^-3trv/7^tr(7~iu;) +2tr(7-^w)^/^)| . (A4) 



This is also linear with respect to the lapse functions. The last one is 



\/ - det g ei{^J g-^ f) = ^/-detg det ^/g-^f = det f = iVdetu . (A5) 



To sum up, the interaction terms read 



$^/3ne„(v^) 

n=0 

N yd^[/3o + f3itT^-f-^u + i/32(trV7"'w - tr(7"iu;)) 



-1, ,\3/2 



-/^af^tr^A/r^- 3trA/r^tr(7"^a;) + 2tr(7"^a;) 
V^det7|/3i + (32tT^fr^ + i/33(^trV7"^w - tr(7~ia;)) | + /3iVdetu 



(A6) 



Notice that all interaction terms are linear in L and A^. This is the advantage of mini- 
superspace model, which makes the Hamailtonian constraint analysis simple. 



In the paper, for simplicity, we always assume that interactions are minimal 2l|, |25| . 
namely 



(3o = 3, (3i = -l, /?2 = 0, f3s = 0, (3^ = 1. 



(AT) 



Clearly, this simplification does not lose any generality concerning with the ghost analysis. 
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Appendix B: Constraint Algebra 



In this appendix, we calculate a Poisson bracket. It is sufficient to look at the following 



VdetcJ ^2 



+ I - ai v/dit^trv/^^, ^= (^P\P^ 3 ' P'^Pij) } 



(Bl) 

The point is that the result is proportional to Oi. Each term can be manipulated as 
T-r- 1 i 7 a 1 ^Vdet 7 9 /I i • 

det 7 , -TT jVr-^,- — TT "'TTj, ^ = — -TT jTT"',- — TT "'TTj, 



^V^det77"'" (7mn7r\ - 27r^^„j 

^^/d^n\ (B2) 



and 



^ / — r- 1 i 7 a 1 9trA/7 9 /I i . ^■ 
tr , -p - p ^p., I = ^ - P ^P, 

= ^p^tr^/^-tr(v/^pa;) , (B3) 

where p represents a matrix with components p"^". 

In this case, Cnl 7^ because there is an interaction between g and /, namely ai ^ 0. 
Thus, if the Poisson bracket is non-trivial or not is determined by the interaction pattern. 
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